Non-local electron transport and cross-resistance peak in NSN heterostructures 
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We develop a microscopic theory describing the peak in the temperature dependence of the non- 
local resistance of three-terminal NSN devices. This peak emerges at sufficiently high temperatures 
as a result of a competition between quasiparticle/charge imbalance and subgap (Andreev) contribu- 
tions to the conductance matrix. Both the height and the shape of this peak demonstrate the power 
law dependence on the superconductor thickness L in contrast to the zero-temperature non-local 
resistance which decays (roughly) exponentially with increasing L. A similar behavior was observed 
in recent experiments. 



Non-local effects in coherent electron transport across 
hybrid structures composed of a superconductor (S) 
attached to two normal terminals (N) have recently 
become a subject of intensive experimentali^i 3 - and 
theoretica l 4 ' i 6 ' 7 ' 8 ' 9 ' 10 ' ' investigations. Provided the 
distance L between two N-terminals (see Fig. 1) is 
smaller than (or comparable with) the superconducting 
coherence length £, two non-local processes contribute to 
electron transport through such NSN devices. One pro- 
cess corresponds to direct electron transfer (DET) be- 
tween two N-metals through a superconductor. Another 
process is the so-called crossed Andreev reflection (CAR): 
An electron penetrating into the superconductor from the 
first N-terminal forms a Cooper pair together with an- 
other electron from the second N-terminal in which case 
a hole goes into the second N-metal. A non-trivial in- 
terplay between DET and CAR yields a rich variety of 
features observed in recent experiments^ "' 




Here we focus our attention only on one of such fea- 
tures, a pronounced peak in the temperature dependence 
of the non-local resistance observed in three-terminal 
NSN structures 1 ^ and attributed to charge imbalance ef- 
fects. Very recently Golubev and one of the authors**- of- 
fered a theory for this phenomenon interpreting the non- 
local resistance peak as a result of a competition between 
the contributions of charge imbalance and Andreev re- 
flection. A striking experimental observation 3 - is that the 
height of the resistance peak depends on the distance L 
between N-terminals much weaker than the correspond- 
ing low temperature cross-resistance which was found 
to decay (approximately) exponentially 1,3 cx exp(— 
in agreement with theoretical prediction o 4 ' 5 ' 6 ' 8 ' 12 . Note 
that due to the restriction L < £ it was not possible to 
address the length dependence of the resistance peaki 3 - 
within the models*. 

Below we will employ the model of three-terminal NSN 
structures with ballistic electrodes 8 - which allows for a 
complete non-perturbative solution of the problem for 
all values of L. We will specifically address the temper- 
ature dependence of the non-local resistance R\2 (T) and 
demonstrate that the height of the charge imbalance peak 



FIG. 1: Schematics of our NSN device. 



scales with L exactly as the corresponding normal state 
resistance Rn 12 . We believe that this observation might 
help to account for recent experimental findings 3 -. In ad- 
dition we will argue that - within the model studied here 
- the charge imbalance peak for JJia(T) can occur only 
in the case of weakly transmitting NS interfaces and it 
quickly disappears as the interface transmissions increase 
beyond the tunneling limit. 

The NSN structure under consideration is depicted in 
Fig. [TJ We will assume that all electrodes are ballis- 
tic and that both NS interfaces (with cross-sections A\ 
and A2) have arbitrary transmissions D\ and D% ranging 
from zero to one. The distance between these interfaces 
L as well as other geometric parameters are assumed to 
be much larger than y/Ai : 2, i-e. effectively both contacts 
are metallic constrictions. At the same time the number 
of conducting channels A/12 = PpAi^/^TV in each contact 
is assumed to be large. 

For convenience, we will set the electric potential of 
the S-electrode equal to zero, V = 0. In the presence of 
bias voltages V\ and V2 applied to two normal electrodes 
(see Fig. [I]) the currents I\ and I2 will flow through 
SNi and SN2 interfaces. These currents can be evaluated 
with the aid of the quasiclassical formalism of noncquilib- 
rium Green- Eilenberger-Keldysh functions^ 3 " g R > A > K . For 
ballistic electrodes considered here the Eilenberger equa- 
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tions take the form 



of our NSN device^: 



era + eV(r, t) - A(r, t), g R ' A ' K (p F ,e, r, t)] + 

+ iv F Vg R - A - K (p F ,e,r,t) = 0, (1) 



where [o, b] = ab — ba, e is the quasiparticle energy, p F = 
mv F is the electron Fermi momentum vector and t 3 is 
the Pauli matrix. The functions g R > A ' K also obey the 
normalization conditions (g R ) 2 — (g A ) 2 = 1 and g R g K + 
g K g A = 0. Here and below the product of matrices is 
defined as time convolution. 

The matrices g R ' A ' K and A are 2x2 matrices in the 
Nambu space 



„R,A,K fR.A.K 
~~ \ fR,A,K ~R,A,K 1 ' 



A 

-A* 



(2) 



and A is the BCS order parameter. Without loss of gen- 
erality below we choose A to be real. The current density 
is related to the Keldysh function g K by the standard 
formula 



j(r,t) 



j de(v F Sp[T 3 g K (p F ,e,r,t)}) 



(3) 



where No — mp F /2tt 2 is the density of state at the Fermi 
level and angular brackets (...) denote averaging over the 
Fermi momentum. 

The above equations should be supplemented by the 
boundary conditions describing electron scattering at NS 
interfaces. Assuming specular reflection at both SNi and 
SN 2 interfaces we introduce their transmission probabil- 
ities Dx, 2 (p Xl , 2 ) = 1 - Ri, 2 (p Xlt2 ) (where p Xl (p X2 ) is the 
component of p F normal to SNi (SN2) interface) and em- 
ploy the standard Zaitsev boundary conditions^ in order 
to match quasiclassical Green functions at both sides of 
each of the two interfaces. Deep inside metallic electrodes 
S, Ni and N 2 the Green functions should approach their 



equilibrium values g 



R,A _ 



±(ef3 — A)/Q, R ' in a super- 



)R,A _ 



conductor and g ' = ±t 3 in normal metals, f2 
v(e± iS) 2 — A 2 . For the Keldysh functions far from in- 
terfaces we have g K = tanh[(e + eVr 3 )/2T](g R - g A ), 
where V = 0, V\ and V 2 respectively in S, Ni and N2 
electrodes. 

The general solution of the problem within the above 
formalism was described in details in Ref.— . Here we 
only point out that the accuracy of the above formalism 
in the case of double-barrier structures under considera- 
tion is justified simultaneously by the two condition o 8 ' 15 : 
A 1>2 « L 2 and jVi >2 > 1. 

At low voltages eVi.2 <C T c we obtain 



h = G n {T)Vi-G X2 {T)V 2 , 
h = -G2i{T)V 1 +G22{T)V 2 . 



(4) 
(5) 



where G\ 2 (T) and GzxiT) are the non-local conductances 



Gia(T) = G 2 i(T) = 



Gn 12 
AT 



de 



cosh 2 (e/2T) 



x (l-^i|a| 2 )(l-^ 2 |a| 2 ) 



1 — tanh 2 iLQ/v F 



(6) 



P(n u n 2 

Here we defined T>\ 2 = 1 — 7£i,2 = ^1 2(f>F7i,2) an< A 
Pf7i(2) is normal to the first (second) interface com- 
ponent of the Fermi momentum for electrons propa- 
gating straight between the interfaces, fl = Vs 2 — A 2 , 
P{Rx,R 2 ) = |1 - R1R2C1 2 - Q[e(l + RiR 2 a 2 ) + Aa(Ri 
R 2 



Q = n^ 1 t&nhiLfl/vp, a = (Q - e)/A, 



Gn 12 — 



%^ 2 MxM 2 V{D 2 
R q p 2 F L 2 



(7) 



is the non-local conductance in the normal state, R q = 
2ir j e 2 is the quantum resistance unit. 

The conductance Gn(T) of the first interface is domi- 
nated by the standard BTK expression^ 

M f (l + |a| 2 )cfe /KJ n l-i?i|a| 2 



Gn(T) 



R q T J cosh 2 (e/2T) \ v F 



while a non-local correction to |5]| is small in the param- 
eter A2/L 2 and will be omitted here. The conductance 
G 22 (T) of the second interface is defined analogously. 

In the temperature interval e~ A / T <C 1 we obtain fol- 
lowing expressions for the conductances 



Gn(T) = Gn(0) + G Nll V2 7 r- l /-e- A / T , 



G 12 (T) 

where 
Gi 2 (0) _ 



2Gjv 12 e 
Gi 2 (0), 



-A/T 



V{D 2 « e~ A / T , 
A 



ln(l/[DiP 2 ])' 



V{D 2 {1 - tanh 2 LA/v F ) 



G Nl2 [1 + TZ 1 Tl 2 + (TZi +K 2 )t&nhLA/v F } 2 '' 



(9) 
(10) 

(11) 



G Nll = 



2Ni 
R„ 



v 



dauer) conductance and Gn(0) 



—^—Di{Pxi) ) is the normal state (Lan- 
v F j 

4A/i / \v. 



Rn 



] -D\{p Xl ) 



v q x V F 

is the subgap (BTK) conductance of the first NS inter- 
face. 

Turning now to the non-local resistance 



Ri 2 (T) = 



G 12 (T) 



(12) 



Gii(T)G 22 (T)-Gi 2 (T)G 2 i(T)' 

we substitute the expressions ([TU]) and © into Eq. (Ql 
and, neglecting the small cross-conductance term G\ 2 G 2 \ 
in the denominator, we obtain 



Ri 2 (T) 



2R Nl2 e- A / T 



Gn(0)/G Nll + y/27rA/Te- A / T 
1 



G 22 (0)/Gat 22 + ^AjTe-^T 



(13) 
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FIG. 2: Non-local resistance R12 of an NSN device as a func- 
tion of temperature for L = vf/T c . Provided both interface 
transmissions are sufficiently low the non-local resistance ex- 
hibits a well pronounced peak. This peak disappears at higher 
barrier transmissions. 



where Rn 12 = Gn 12 / '(Gjv^Gjv 22 ) i s the cross-resistance 
in the normal state. Eq. (|13p applies within the interval 
-D1-D2 <C e~ A / T -C 1 and represents the central result of 
this paper. 

We observe that the non-local resistance (fl~3)) first 
increases with decreasing T reaching its maximum at 
T = T* and then decreases with T at lower temperatures. 
For the temperature T* with the logarithmic accuracy we 
obtain 



hi 



Gnxx^N22 

Gn(0)G 22 (0) 



(14) 



or simply T* ~ A/ln(l/£>) for symmetric structures 
with D1.2 = D. This result matches qualitatively with 
that obtained in Ref—^ within a different model. 

The temperature dependence of the non-local resis- 
tance i?i 2 (T) for our NSN device is depicted in Fig. 
2 for different values of the interface transmissions. In 
the limit of low transmissions the resistance i?i 2 (T) de- 
creases with T right below the superconducting critical 
temperature T c but then turns upwards and exhibits a 
well pronounced peak. At lower T the resistance Rti(T) 
decreases sharply and eventually tends to i?i 2 (0) = 
Gi 2 (0)/[Gn(0)G 22 (0)] in the limit T -> 0. With in- 
creasing interface transmissions £>i )2 the peak gets less 
pronounced and eventually disappears, in which case the 
non-local resistance i?i 2 (T) decreases monotonously with 
temperature. The same feature can also be observed in 
Eq. (HU). 

As it is clear from Fig. 3 both the height and the form 
of the non-local resistance peak scale with the distance 
L exactly as the normal state value Rn 12 - For instance, 



Dx = 0.1, D 2 = 0.1, L = vf/T c 
Dx = 0.03, D 2 = 0.03, L = vp/T,, - 
Dx = 0.03, D 2 = 0.03, L = 5vf/T c - 
Dx = 0.03, Di = 0.03, L = Wv F /T c ■ 




FIG. 3: The same as in Fig. 2 for different values of the 
distance L between two N-terminals (measured in units of £). 
Three curves corresponding to the same barrier transmissions 
but different values of L practically coincide, i.e. in the vicin- 
ity of the peak R\2(T) in the superconducting state scales 
with L exactly as Rn 12 - 



from Eq. (fT3"|) we easily determine the height of the re- 
sistance peak i?i 2 (T*) which reads 



Rl2(T* 

RNxi 



2T* 

7rA 



/Gn(0) ^ 


/G 22 (0) 


/ Gnxx 


Y Gat 22 



(15) 



We observe that the right-hand side of Eq. (|15|) is ex- 
pressed only in terms of local conductances and, hence, 
is independent of L. This is an important result which 
might account for experimental observations^ of a much 
weaker L-dependence of i?i 2 (T*) as compared to the 
zero temperature value i?i 2 (0). For instance, within our 
model of ballistic electrodes for L > (0 we have i?i 2 (0) tx 
cxp(-2LA(0)/wf) while R 12 {T*) cx 1/L 2 . Even weaker 
length dependence of i?i 2 (T*) oc i?Ar 12 oc 1/L is expected 
in the diffusive limit. It would be interesting to scale the 
data- for the resistance peak at different lengths L with 
the corresponding normal state resistance Rn 12 (L) ■ 

Although our theory correctly accounts for some key 
features of the experimental data it is important to bear 
in mind that the model employed here deals with ballistic 
electrodes connected via metallic constrictions whereas in 
experiments^ the electrodes were most likely diffusive. 
In addition, we disregarded any relaxation mechanisms 
for non-equilibrium quasiparticles inside the supercon- 
ductor (except for their escape into the normal terminals) 
while such mechanisms (caused, e.g., by electron-phonon 
and electron-electron interactions) are obviously present 
in experiments being responsible for a finite charge im- 



4 



balance length. Hence, one can also expect certain differ- 
ences. For instance, within our model the non-local resis- 
tance peak occurs only at small interface transmissions 
while the authors^ observed such a peak at moderately 
high transmissions of NS interfaces. This might be an in- 
dication to relatively more pronounced charge imbalance 
effects. 

In summary, we have developed a quantitative the- 
ory describing the peak in the temperature dependence 
of the non-local resistance recently observed in three- 
terminal NSN devices^ at sufficiently high T. This 
peak emerges as a result of a trade-off between quasi- 



particle/charge imbalance and subgap (Andreev) contri- 
butions to local and non-local conductances of the device. 
Both the height and the shape of the peak scale with the 
normal state resistance Rn 12 thus demonstrating much 
weaker (power law) dependence on the superconductor 
thickness L as compared to the zero-temperature resis- 
tance i?i2(0) which decays (roughly) exponentially with 
increasing L. 
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